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Abstract 

Propagation of electromagnetic plane waves in some directions in gravitationally 
affected vacuum over limited ranges of spacetime can be such that the phase veloc- 
ity vector casts a negative projection on the time-averaged Poynting vector. This 
conclusion suggests, inter alia, gravitationally assisted negative refraction by vacuum. 

Keywords: General theory of relativity, Gravitational lensing, Hidden matter, Negative 
phase velocity, Poynting vector, Special theory of relativity 

1 Introduction 

The discovery of (purportedly) isotropic, homogeneous, dielectric-magnetic materials that 
bend electromagnetic rays the "wrong way" [1] created quite a stir in 2001 [2], with claims 
and counterclaims flying all around [3, 4]. The situation has recently been settled, with 
unequivocal demonstrations by several independent groups [5]-[7]. See Ref. 8 for a compre- 
hensive review. A range of exotic and potentially useful phenomenons — such as negative 
refraction, negative Doppler shift and inverse Cerenkov radiation — have been predicted for 
materials of this type, wherein the phase velocity is directed in opposition to the energy 
velocity as quantified through the time-averaged Poynting vector. These materials have sev- 
eral names, including left-handed materials, negative-index materials, and negative-phase- 
velocity (NPV) materials. We prefer the last term [8]. 

Subsequently, the possibility of NPV propagation of light and other electromagnetic 
waves was established in a variety of anisotropic materials [9]— [11]. In these materials, 
NPV propagation is indicated by the projection of the phase velocity on the time-averaged 
Poynting vector being negative. 
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Even more interestingly, materials that do not permit the observation of NPV propaga- 
tion by observers in a relatively stationary (i.e., co-moving) inertial reference frame have 
been shown, after the invocation of the postulates of special theory of relativity (STR), to 
allow observation of NPV propagation in other inertial frames [12]. That permits one to 
envisage STR negative refraction being exploited in astronomical scenarios such as, for ex- 
ample, in the remote sensing of planetary and asteroidal surfaces from space stations. Quite 
possibly, space telemetry technologies will be the first to reap the benefits of STR nega- 
tive refraction. Application to remotely guided, extraterrestrial mining and manufacturing 
industries can also be envisioned. Furthermore, many unusual astronomical phenomenons 
would be discovered and/or explained via STR negative refraction to interpret data collected 
via telescopes [13]. 

As is well known, vacuum (i.e., matter-free space) appears the same to all inertial ob- 
servers [14]. Therefore, as a co-moving observer cannot deduce the occurrence of NPV 
propagation in vacuum, neither can any observer moving with a constant velocity. This 
could lead one to believe that NPV propagation is impossible in huge expanses of interstel- 
lar space. However, gravitational fields from nearby massive objects will certainly distort 
electromagnetic propagation, which is a principal consequence of the general theory of rel- 
ativity and is indeed used nowadays in GPS systems. Our objective here is to establish 
that gravitationally affected vacuum can support NPV propagation, at least in spacetime 
manifolds of limited extent. 



A gravitational field curves spacetime, which effect is captured through a metric g a p- Elec- 
tromagnetic propagation in gravitationally affected vacuum may be described in terms of 
propagation in an instantaneously responding medium in flat spacetime [15, 16], at least in 
spacetime manifolds of limited extent. That is, the nonuniform metric g a p can be locally ap- 
proximated by the uniform metric g a p [17]. On asssuming the convention g a p = (+, — , — , — ), 
the constitutive relations of vacuum in the equivalent flat spacetime are expressed in Gaus- 
sian units as [16] 
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where ee mn is the Levi-Civita tensor, and 
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(3) 




(4) 



with g = det [<?«,«]. We note that the metric g a p is real symmetric [18]. 



jRoman indexes take the values 1, 2 and 3; while Greek indexes take the values 0, 1, 2, and 3. 



2 



The constitutive relations (1) and (2) can be expressed in 3-vector form as 



D = e 2' E-Tx H, (5) 
B = M-I + ExE, (6) 

wherein SI units are implemented. The scalar constants e and /i denote the permittivity and 
permeability of vacuum in the absence of a gravitational field, respectively, and c = ^/l/e /x . 
Our coordinate system is chosen such that the second-rank cartesian tensor 7 is diagonal; i.e., 

7 = diag (7 X , 7 y , 7 2 ). In the gravitational field of a mass rotating with angular momentum J, 
the gyrotropic vector T is proportional to R x J, where i? is the radial vector from the centre 
of mass to the point of observation. By considering a small region of space at a sufficiently 
remote location from the centre of mass, we take £ to be independent of R. 
We seek planewave solutions 

E = Re {E exp [i (fc • r - ut)] } , (7) 
# = Re{H exp[i(k -r-cut)]} , (8) 

to the source-free Maxwell curl postulates 

Vx£ + ^B=0, (9) 

Vx /J -SLd = 0. (10) 

Here k is the wavevector, r is the position vector, u is the angular frequency, and t denotes 
the time; whereas E and H_ are complex-valued amplitudes. 

An eigenvector equation for E_ is developed as follows. By combining (5)-(8) with the 
Maxwell curl postulates, we derive 

pxE = ujfi g-H , (11) 
pxH = -ue g'E , (12) 

in terms of 

p = k-uT. (13) 
The use of (11) to eliminate i? from (12) provides, after some simplification, 

K'E o = 0, (14) 

where 

V\^= (^k^det 7 — p"y-p^jl + pp"y, (15) 

and the notation k = u^e /i has been introduced. A dispersion relation thus emerges from 
(14) as 

det [m] =0, (16) 
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which may be expressed in the form 



k det 



7 



k det 



7 



0. 



Hence, we conclude that planewave solutions satisfy the condition 



(17) 



(18) 



p • 7 • p — kl det 7 

Let us consider eigenvector solutions to (14). Substitution of (18) into (14) provides 

PP'1'E = 0; (19) 

thereby, all eigenvector solutions are necessarily orthogonal to p • 7. To proceed further, let 
us — without any loss of generality — choose the wavevector k to lie along z axis, and the 
vector T to lie in the xz plane; i.e., 



k — kzt z , 

£ = r (u x sin 9 + u z cos 9) , 



(20) 
(21) 



where u x , u y and u z are the cartesian unit vectors in the equivalent flat spacetime. Since 

p*7 = — ouTj x sin 9u x + (k — 00T cos 9) 7 2 m 2 , (22) 

it is clear that two linearly independent eigenvectors satisfying (18) may be stated as 

ei = %, (23) 
e 2 = u y x (p • 7) (24) 
= (k — uT cos 9) ^zikx + sin 6*m 2 . (25) 

After assuming that 7 is invertible, we deduce the corresponding magnetic field eigenvectors 
from (11) as 

1 



h 2 = 



LVH = 

1 



7 1 • [(cuF cos 9 — k)u x — ujT sin ] , 



Hence, the general solution is given by 



7 1 * \{k — ujY cos9) 7 Z + (c^r sin 9) 7 X ] 



+ C 2 [(A; - c^r cos 9) ^ z u x + cjr7 x sin 9u z ] , 
' Ci [(c^r cos 9 — k) u x — ujT sin 9u z ] 



(26) 
(27) 

(28) 



-1 



ua - 



7 



+C 2 [(k - ujT cos 9f lz + (c^r sin 9) 2 lx ] u y \ , 



(29) 



wherein C\ and C 2 are arbitrary constants. 

The wavenumbers arise from the dispersion relation (17) as follows. Substituting (13) 
into (18), we obtain the /c-quadratic expression 



k 2 ^ z - 2k^ z uY cos 9 + uj 2 T 2 sin 2 9 + ^ z cos 2 9) - k 2 det 
since k • £ = kT cos 9. The two /c-roots of (30) are 

k + = uj [y cos 9 + ^Je fj, j x j y — —T 2 sin 2 9^j , 
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k' 



lu (v cos 9 - y eoHo'jx'jy - ^r 2 sin 2 6^ 



Finally, let us consider the time-averaged Poynting vector given by 

P = \Re{E xH* }. 



(30) 

(31) 
(32) 

(33) 



After utilizing the general solution (28) and (29), the component of the Poynting vector 
aligned with the u z axis is obtained as 

u z -P = — — {k-ujTcos9) (\d\ 2 + \C 2 \ 2 "f z uj 2 det \-f]) . (34) 

2w/i 7 2 ^ L = J/ 

The energy density flow in the direction of the wavevector k + , corresponding to the root 
k + given in (31), is thus 

1 



P = 



7 



7, 



2/i,,7. 

x (\d\ 2 + \C 2 \ 2 lz u 2 det [2]) 



r cos 9 J eo/ioTaTj, - — r 2 sin 2 6» + ( e /i 7s7i/ - — r 2 sin 2 9 



lz 



Let us notice that the inequality 



toUolxly r sin 9 > 

7* 



(35) 



(36) 



must be fulfilled in order for k to be real-valued. Therefore, the defining inequality for NPV 
propagation, namely 

k + ' P < , (37) 

is satisfied provided that 



r cos 9 > x /e ^o7a!7i/ ~ — r2 sin2 

lz 



(38) 



holds. Analogously, we find that NPV propagation is signalled for the k wavenumber by 
the condition 



r cos 9 > J eo^Q^y - — r 2 sin 2 9. 
V 7z 



(39) 



In deriving (38) and (39), we used the fact that l x ,y,z > by virtue of the signature of g a p. 
We note that the NPV conditions (38) and (39) are independent of frequency. 
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3 Concluding Remarks 



The inequalities (38) and (39) can be satisfied for specific ranges of the angle 9, for given 
7 and T. Thus, we have shown that NPV propagation in some directions is possible in 
gravitationally affected vacuum over limited ranges of spacetime. The possible existence of 
gravitational fields which can deliver 7 and T necessary for the satisfaction of (38) and/or 

(39) is a matter for astrophysicists to discuss. 

We are content here to state that, just as scientific and technological applications of 
STR negative refraction (by materials) can be envisaged [12, 13], similar and different conse- 
quences of gravitationally assisted negative refraction by vacuum are possible. In particular, 
designers of channels for space communication shall also have to account for the possibility 
of negative refraction due to massive objects between the two ends of every channel. 

Furthermore, current ideas on the distribution of mass in the as-observed universe may re- 
quire significant revision, since electromagnetic signals from distant objects may be deflected 
in a manner which has not hitherto been accounted for. Thus, our work has implications for 
gravitational lenses [19]. Gravitational lensing involves nonuniform metrics, and the distri- 
bution of matter in the universe has been constantly changing. While the spatiotemporally 
local evolution of the universe may be deduced adequately from electromagnetic signals re- 
ceived by our telescopes, reasonably accurate deductions about the spatiotemporally global 
evolution of the universe from similar measurements may be particularly difficult to make - 
owing to gravitationally assisted negative refraction. 

Finally, our work suggests that it is time that research on consequences of gravitationally 
assisted negative refraction by materials be undertaken. 
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